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Superfluid — Bose glass transition in two dimensions 
at T = 0: analytic solution of a mode-coupling toy model. 

E.V. Zenkoyi'B 

^ Ural State Technical University, 620002 Ekaterinburg, Russia 

Analytic expression for the memory function and the optical conductivity of the two-dimensional 
Bose gas with logarithmic interaction at T = in presence of point-like impurities is obtained within 
the mode-coupling approximation. Depending on the value of a dimensionless combination of the 
model parameters proportional to the strength of the impurity potential, two different phases are 
distinguished, viz. the disordered superfluid and insulator (Bose glass), separated by an intermediate 
(quasi)metal state. 

Keywords: glass transition, memory function, mode-coupling theory, Bose condensate 



I. INTRODUCTION 

The behaviour of many-particle classical and quantum 
systems in presence of a random perturbation is a long- 
standing problem of great scientific interest. 

Various methods have been developed to describe the 
behaviour of a system subject to the random external 
force, such as the coherent potential and other effective 
medium types of approximationsii^, the self-consistent 
hydrodynamic scheme for the description of the over- 
damped charge dynamics in presence of a quenched disor- 
der, elaborated by W. Gotze and coworkers^, the replica 
theories of flux lattices and Wigner crystals^. 

It is interesting to observe, that for very different model 
systems in the localization regime, many quantities of 
physical interest, such as the dynamic conductivity, were 
found to manifest an intimate resemblance. E.g., in the 
three-dimensional (3D) electron system, the Drude tail 
of the optical conductivity flattens out and gradually 
evolves into a flnite energy peak as the system approaches 
the Anderson transition^ from the metallic side of the 
phase diagram. Similar results are obtained for the two- 
dimensional (2D) electron gas and also for 2D and 3D 
Bose gas at zero temperature^iZ, where the effects of the 
disorder lead to the decay of the coherent response of the 
condensate acoh and the transfer of the spectral weight 
to a finite-energy broad feature. The same general pic- 
ture holds also for pinned density waves^ and Wigner 
crystals^, where the conductivity displays the inhomoge- 
neously broadened feature at the pinning energy. 

However, the understanding of the disorder effects in 
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complex quantum systems is still an issue, that have been 
given the new impetus since the discovery of the strongly 
correlated oxides. An extreme complexity of these prob- 
lems usually rules out the possibility of a tractable an- 
alytic treatment. In this article we present a simple ex- 
actly solvable model of the superfluid-insulator transition 
in the 2D Bose gas with logarithmic interaction in pres- 
ence of the random short-range potential of the point 
defects at zero temperature. 

The key result of our analysis is the analytic expression 
for the memory function, obtained within the framework 
of the mode-coupling theoryS. This enables us to ex- 
plore the evolution of the optical conductivity in the full 
range of the relevant physical parameters and to trace the 
superfluid-insulator transition, that takes place, when a 
dimensionless parameter of the disorder exceeds a critical 
value. Despite a number of simplifying assumptions, we 
believe, that the theory, developed in the next section, 
is of a certain interest, being an elementary but instruc- 
tive exactly solvable model of the disorder-induced glass 
transition. 



II. THE MODEL 

The central role in our analysis plays the memory func- 
tion M{z) = M'+iM", that measures the relaxation rate 
of the current excitations and incorporates the effects of 
an external random perturbation due to the impurities. 

Within the mode coupling theoryi the memory func- 
tion obeys the self-consistency relation: 



M{z) = i'y + 



%ik,z + M{z)) 



2nm J 1 + M(z)$o(fc, z + M(z))/x(fc, 2 = 0) 



fc2(|i7(fc)p)dk. 



(1) 



where 7 is a bare relaxation constant, corresponding to 
some internal current decay channels, U{k) is the Fourier 
transform of the random potential, n and m are the con- 



centration and the mass of the particles, respectively. 
The angular brackets stand for the averaging over the 
disorder. The density response function $0 is related to 
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the generalized susceptibility xo{k, z) 
out disorder: 



of the system with- 



Xo{k,z) - xo{k,0) 



(2) 



Thus, the free model being specified by its correlation 
functions, eq. (H]) offers a direct way to explore the effects 
of the disorder on its electric and optical properties. 

We consider the 2D Bose gas at zero temperature, 
where all particles are condensed into the superfluid 
phase. Hence, no bare mechanisms of relaxation are 
presumed2£ and 7 = in ([1]). The susceptibility of non- 
interacting 2D bosons is readily obtainedii: 



Xnonint . ^) 



2nek 



(3) 



where Sk = k'^ /2m is the dispersion of free particles. The 
interparticle interaction V modifies the susceptibility. A 
qualitative estimate of this effect may be done within the 
random phase approximation: 



Xo{k,z) = 



X.nonint. 

{k,z) 



1 + V{k)Xnonint.ik,z)' 



(4) 



Having in mind to work out a simple tractable model, we 
neglect here a certain limitation of this mean-field type of 
approximation in low dimensions, as well as the local-field 
corrections^^'. On the other hand, the mean-field concept 
retains its significance as a highly intuitive approach, that 
makes it possible to get insight even into the problems, 
where its formal substantiation is problematic, e.g. in the 
physics of the superfluid filmsi^ and the plane vortices in 
quasi-two-dimensional magnetsiW, 

In purely two-dimensional systems, such as the planar 
vortices, the axially symmetric solution of the 2D Laplace 
equation provides the following form of the interparticle 
interaction: 



V{r)^QHn[^) 



(5) 



(the two-dimensional analog of the Coulomb potential), 
where Q is the boson charge and tq is an appropriate 
length scale, that may be considered as a mean interpar- 
ticle distance, vq ~ n^^^, where n is the two-dimensional 
concentration of particles. Thus, the potential V{r) is re- 
pulsive for r < tq and attractive for r > tq, favouring the 
homogenisation of the system. The ground-state proper- 
ties of the free 2D Bose gas with the interaction ([5]) have 
been investigated Correspondingly, the potential in 
momentum space is: 



V{k) 



2 7rQ2 
k^ 



(6) 



that may be obtained as the 2D Fourier transform of 
V(r) exp(— a k r) in the limit a — ^ 0. 

Hereafter it is convenient to introduce^he dimen- 
sionless variables p = p/po, lo = ^l^o, M = M/eq, 



defining the units of momentum and energy as po = 
(STrQ^ft^mrt)-^/'* and Eq — Po/^™; respectively^!. In what 
follows, we shall consider the case of identical impuri- 
ties, with random coordinate distribution. The average 
((. . . )) in rhs. of cq. ([T]) then reduces to the multiplica- 
tion by the impurity concentration rii. Substituting eqs. 

with ^ into ([T]) yields the nonlinear equation on 
the memory function, where all quantities are reduced to 
the dimensionless form and tilda's are omitted: 



oc 

f (z + M)fcV(fc)2 

TT J (l + fc4)(l-|-fc4-z(z + M)) 




dk. (7) 



Thus, the optical response of the model is governed by a 
single parameter 



A - -— I — I - "'"^ 
2 n \eo J 



El 

HQ 



(8) 



where Uq measures the strength of the impurity poten- 
tial, U{k) = Uo(p{k). Note, that Uo has the dimension of 
energy, and that of </'{fc) is The similar non-linear 

equations on M{z) have been derived early on for the 2D 
and 3D Bose gas with the Coulomb repulsion^'^. How- 
ever, their analytic solution were never possible because 
of the complexity of the integrand. 

In what follows, we focus on the special case of the 
short-range impurity centers, that is usually considered 
in the theory of pinned density waves and Wigner crys- 
tals: f{r) — S(r), whence f{k) = 1. Then, the integral in 
rhs. of ([7]) can be solved to yield the quadratic equation 
for M , whence we obtain the explicit expression for the 
memory function: 



M = 2A I 



1 - A 



•f - A 



- 1 



(9) 



where the root with the non-negative imaginary part is to 
be chosen. This is the main result of the present article. 
Note, that within the mode-coupling theory eq. jl]) is an 
exact relation. 



III. OPTICAL CONDUCTIVITY AND THE 
SUPERFLUID - GLASS TRANSITION 

The optical response of the system is most conveniently 
studied in terms of the optical conductivity. Within the 
memory function approach this is expressed in the gen- 
eralized Drude form: 



'z + M(z)' 



(10) 



where z = lu + Oi is the complex frequency, cOp is the 
plasma frequency (w^ = Q^n/m, where n, Q and m are 
the concentration, charge and mass of bosons, respec- 
tively). 
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FIG. 1: The optical conductivity spectra of 2D Bose gas with 
point- like impurities at different levels of the disorder: 1 - A = 
0.2 (superfluid), 2 - A = 1.0 (metal), 3 - A = 1.8 (insulator). 
The arrow indicates the edge of the incoherent wing, that falls 
to the same frequency for both curves 1 and 3. For the sake 
of clearness, the infinitely narrow 5-like peak in the superfluid 
phase is shown broadened. 



The explicit analytic expression for the memory func- 
tion ^ makes it possible to explore various regimes of 
the conductivity of the system. First of all, it is straight- 
forward to show, that this results preserves the sum-rule: 



a'{uj)duj = ^TTWp, 



(11) 



that we notice as an important sign of the sound be- 
haviour of the present model. 

Turning to the detailed analysis of the conductivity 
spectrum, let us consider first the weak disorder limit, 
A <C 1. Expanding the generalized Drude formula (fTO|) 
to first order in A, we arrive at the following expression 
for the real part of the optical conductivity at finite fre- 
quencies: 



C''{uj)incoh — 



2A 

0, 



1 



UJ <l. 



(12) 



Thus, the spectrum shows a gap, that a more accurate 
treatment finds to extend at = 1 1— A| . The free system 
parameters being fixed, this gap frequency is governed by 
the amplitude of the disorder (pinning) potential. Above 
the gap, the optical conductivity cr[ncoh displays a broad 
asymmetric feature, shown in fig. |TT1 It corresponds to 



S(coo) 




FIG. 2: The dependence of the spectral weight S of the inco- 
herent wing, ea. (|15|l . on the disorder parameter A. The knee 
at Ac = 1 indicates the superfluid-insulator transition. 



the localized plasma oscillations of the charge density 
and is hereafter called the incoherent wing. Such a be- 
haviour is typical of the charge density waves. Interest- 
ingly, we recover in p2p the expression for the conduc- 
tivity of pinned CDW, first derived by Fukuyama and 
LeaiS,. We believe this result to be the consequence of 
a nontrivial analogy between the pinned CDW and the 
disordered Bose condensate, albeit such a literal coinci- 
dence of the formulas is rather accidental. The point is 
that pinned CDW is a collective state, that can be de- 
scribed by a model of a particle in a potential welU^. On 
the other hand, the condensate at T = is also a collec- 
tive (coherent) state, that behaves as a single particle. 

However, it follows from the sum- rule pTjl . that the 
contribution of (Tmcoh (US) to the optical conductivity 
spectrum cannot be the only one, since it vanishes as 
A tends to zero. As ea.([T^ is reasonable at w large, the 
careful examination of the conductivity at low frequencies 
is required. The corresponding expansion of the memory 
function (O is: 



M{z) ~ {\ - sign{\ - K)) 



A(l-A) 



A 



1 - A 



sign{\- 



-A)z. 
(13) 

The behaviour of the system depends critically on the 
magnitude of A. According to eq. p3p . the memory 
function is linear near the zero frequency for small dis- 
order parameter, (A < 1). The corresponding optical 
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conductivity is: 



Cr{LL!)coh = (1 - A)- 



2— OJ+Oz 



(14) 



This expression describes the coherent response of the 
condensate. Its spectral weight drops hnearly as the dis- 
order parameter A increases. According to the sum-rule 
pTjl , this implies the emergence of the new spectral fea- 
ture at higher frequencies, described by Gmcoh (US)- The 
full conductivity spectrum is the sum of the coherent (fH|) 
and incoherent (|12p contributions, as depicted in fig. |TT1 

In the strong disorder regime (A > 1), one obtains 
from in the leading order M ^ — 1/z, and (t'(0) = 
0. Within the theory of the liquid-glass transitions^^ 
this pole of the memory function is known as the non- 
ergodicity pole and its appearance is identified with the 
onset of the glass phase as it calls forth the long-lived 
density fluctuations {Sp{k,t)Sp{k,0)) , that do not decay 
to zero at large times. Thus, we call the insulating phase 
of our model the Bose glass. 

The superfluid-glass transition may be clearly visual- 
ized by considering the partial spectral weight of the in- 
coherent contribution. 



^incoh 



a'{uj)duj. 



(15) 



1-A| 



This quantity is shown in figH] against the disorder pa- 
rameter A. It can be seen, that S grows progressively as 
A approaches unity from below, i.e. the spectral weight 
moves away from the condensate to the localized modes. 
However, it does no longer change as soon as A > 1, al- 
though ujQ continues to grow. This is because no spectral 
weight is now contained below uq and the superfluidity 
is completely suppressed by disorder. 

The system undergoes the superfluid-insulator transi- 
tion via an intermediate state, that may be called metal- 
lic. The high-frequency asymptotics of the optical con- 
ductivity in both phases follows the Drude law with the 
damping parameter equal to 2A, as can be seen from 
the limiting form of At the critical transition point, 
A = 1, the spectral weight of the condensate vanishes, 
and so does the gap frequency wq, while the memory 
function is M = 2z. This obviously corresponds to the 
conventional Drude optical conductivity with the disper- 
sionless relaxation rate in the full spectral range. 



IV. SUMMARY AND CONCLUSIONS 

In this paper we have studied analytically the 
superfluid-glass transition within the purely two- 
dimensional Bose gas with logarithmic interparticle in- 
teraction at T = K with static disorder, represented by 
the point-like random impurities. The dielectric proper- 
ties of the system are considered proceeding from the ex- 
plicit expression for the memory function M{z), obtained 
within the mode-coupling theory. The onset of the glass 



phase manifests itself as the appearance of the nonergod- 
icity pole^^ of the memory function, M(z) ~ 1/z. 

The above analysis shows, that the phase diagram of 
the model is delimited by the following combination of 
the physical parameters: 



4n2 yhQ 



(16) 



Untill A < 1, the optical conductivity cr(w) spectrum 
consists of the response of the condensate oc S{uj), where 
S{w) is the Dirac delta function, and an inhomogeneously 
broadened feature at w > |1— A|. The latter one is typical 
to the pinned density waves and falls off as 1/ cj^ at large 
frequencies. 

At A > 1 the Dirac delta function contribution disap- 
pears and the system undergoes the transition to a local- 
ized phase. Hence, large value of either the boson effec- 
tive mass m, the amplitude of the impurity potential Uq 
or the concentration of impurities n^, other parameters 
being fixed, favours the localization and suppresses the 
superfluid phase, while strong boson interparticle repul- 
sion Q and large boson concentration n have an opposite 
effect. 

A natural question arises on the internal structure of 
the novel localized phase. Within the present hydro- 
dynamic approach this problem cannot be resolved on 
the microscopic level. Qualitatively, the phase transition 
may be realized as a condensate fragmentation in the 
rugged potential landscape, that entails the breakdown 
of the long range non-diagonal order. The system breaks 
up into the patches, trapped by the randomly distributed 
pinning centers. By analogy with the mode-coupling the- 
ory of simple liquids, where the similar transition is re- 
ferred to as the liquid-glass transition, we call the local- 
ized phase of the present model the Bose-glass phase (see 
alsc^'^2, where this name was first introduced). 

Now, let us dwell in more details on the comparison 
of our results with some other approaches to the local- 
ization and superfluid-insulator transition in a system of 
2D bosons. In Rcf.^ the phase diagram of the disordered 
2D Bose liquid is discussed in terms of two dimensionless 
parameters, J/V and A/l^, where J is the zero point en- 
ergy of the bosons, V measures their short-range repul- 
sion and A is a disorder parameter. It was found, that 
large value of J/V favours the superfluid (SF) phase, 
while the decreas of J/V or the increase of A/V leads 
to the localization and the supression of SF phase. In 
our case the only parameter A, eq.®, governs the phase 
diagram. The ratio n^/rii, that enters therein may be 
understood as the concentration of bosons per a patch 
around an impurity. Thus, the first factor in A is propor- 
tional to the inverse energy per boson in such a patch^-, 
niiri/ {h?n'^) ~ J^^. Hence, A may be written in terms 
of Ref.-2. as: A ~ A'^/{JV^). It can be seen, that the 
increase of A, that describes the gradual supression of 
SF phase and the onset of the localization, qualitatively 
agrees with the mentioned general results of Ref;^. 
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The problem of the relevance of the model to the re- 
alistic physical systems is subtle, since the real parti- 
cles, either boson or fermions, are likely to interact via 
not logarithmic, but rather the screened Coulomb po- 
tential. However, the effective model of 2D Bose liquid 
with the logarithmic interaction at T = is related to 
the physics of the pinned vortices in the bulk type II 
superconductorsi^,. We do not scrutinize this question 
any further in the present article. 

We believe the model considered to be of importance 
on its own right, as it exhibits a rather non-trivial be- 
haviour and allows a fully analytic treatment. The 
comparison of our results with those, obtained for the 
other model of the disordered Bose systems numeri- 
cally and semi-analytically within the same theoretical 
framework^, leads to the conclusion, that in both cases 



the phase diagram and the optical spectra are essentially 
the same and are fairly weakly sensitive to the charac- 
ter of the interparticle interaction, the impurity potential 
and the dimensionality of the systems. Thus, we hope the 
present study may shed light on the behaviour of more 
realistic models, that do not permit such an elementary 
analysis. 
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